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Abstract 

In the centennial of Ettore Majorana's birth (1906 — 1938?), we 
re-examine some aspects of his fundamental scientific production in 
atomic and molecular physics, including a not well known short com- 
munication. There, Majorana critically discusses Fermi's solution of 
the celebrated Thomas-Fermi equation for electron screening in atoms 
and positive ions. We argue that some of Majorana's seminal contri- 
butions in molecular physics already prelude to the idea of exchange 
interactions (or Heisenberg-Majorana forces) in his later works on the- 
oretical nuclear physics. In all his papers, he tended to emphasize the 
symmetries at the basis of a physical problem, as well as the limita- 
tions, rather than the advantages, of the approximations of the method 
employed. 

Key words: Ettore Majorana; Enrico Fermi; Thomas-Fermi model; 
exchange interactions; atomic and molecular models; neutron. 



1 Introduction 

Ettore Majorana's most famous, seminal contributions are certainly those 
on the relativistic theory of a particle with an arbitrary instrinsic angular 
momentum pQ, on nuclear theory j2], and on the symmetric theory of the 
electron and the positron [3] . In particular, the latter paper already contains 
the idea of the so-called Majorana neutrino as has been correctly em- 
phasized [4 . The quest for Majorana neutrinos is still the object of current 
fundamental research (see, e.g., Ref. [S], and Ref. [S] for a general overview). 

In this note, we would like to reconsider two more papers by Majorana 
IE] , both on atomic and molecular physics, and show how they are pre- 
cursor to his theoretical work on the exchange nuclear forces, the so-called 
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Heisenberg-Majorana forces [2|. We will also try and emphasize his criti- 
cal sense and great ability to catch the relevant physical aspects of a given 
problem, beyond his celebrated mathematical skills, as witnessed by contem- 
poraries and colleagues who met him personally Ell HI] (see especially 
Ref. for more references). 

Both Amaldi and Segre have provided us with a vivid account of Ma- 
jorana's first meeting with Enrico Fermi. Majorana and Fermi first met in 
1928 at the Physical Institute in Via Panisperna, Rome. At that time, Fermi 
was working on his statistical model of the atom, known nowadays as the 
Thomas-Fermi model, after the names of the two authors who derived it in- 
dependently [131 CH CHI • Such a model provides an approximate alternative 
to solving Schrodinger equation [TH] , and paved the way to density functional 
theory [T7] . 



2 Thomas- Fermi model 

Within Thomas-Fermi approximation, the electronic cloud surrounding an 
atom is described in terms of a completely degenerate Fermi gas. Following 
Ref. JHl, one arrives at a local relation between the electron density p(r) 
at position r with respect to the nucleus, and the momentum pf(t) of the 
fastest electron (Fermi momentum), as 

Air 

p(r) = 2- t 4(t), (1) 

where the factor of two takes into account for Pauli exclusion. In Eq. 
the Fermi momentum pp( r ) depends on position r through the self-consistent 
potential V(r) as 

p 2 F (r) = 2m[E F -V(r)}, (2) 

where E F is the Fermi energy, and m is the electron mass. Fermi energy E-p 
is then determined via the normalization condition 

J p(r) rf 3 r = N, (3) 

where N is the total electron number, equalling the atomic number Z for a 
neutral atom. Inserting Eq. (J2J) into Eq. ((Tj), making use of Poisson equation, 
and introducing Thomas-Fermi screening factor through 

Ze 2 

V(r)-E F = 0(r), (4) 
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one derives the adimensional Thomas- Fermi equation for a spherically sym- 
metric electron distribution, 



where 



and b sets the length scale as 



dx 2 XV2' ^ 



r = bx, (6) 



1 /9tt 2 \ 1/3 0.8853 



b = 4 \JZ) ao = ^7^ G °' (?) 

with a the Bohr radius. 

Eq. (jSJ) is 'universal', in the sense that the sole dependence on the atomic 
number Z comes through Eq. ((7|) for b. Once Eq. (jSJ) is solved, the self- 
consistent potential for the particular atom under consideration is simply 
obtained by scaling all distances with b. 



3 Asymptotic behaviour of the solution to 
Thomas-Fermi equation 

Fermi endeavoured to solve Eq. (j5J) analytically without success. On the 
occasion of his first meeting with Majorana, Enrico succintly exposed his 
model to Ettore, and Majorana got a glimpse of the numerical results he had 
obtained over a week time, with the help of a primitive calculator. The day 
after Majorana reappeared and handled a short note to Fermi, where he had 
jotted down his results. Majorana was amazed that Fermi's results coincided 
with his own. 

How could Majorana solve Eq. (jSJ) numerically in such a short time with- 
out the help of any calculator? Various hypotheses have been proposed. Did 
he find an analytical solution? At any rate, there are no physically accept- 
able analytical solutions to Eq. (jSJ) in the whole range < x < +oo. The 
only analytical solution, 

144 

= -jT> ( 8 ) 

would have been found later by Sommerfeld in 1932 ^H], and is physically 
meaningful only asymptotically, for x 3> 1. 
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The most likely hypothesis is probably that of Esposito (THj, who, together 
with other authors 20J, has found an extremely original solution to Eq. (J3J) in 
Majorana's own notes (see also Ref. [21])- The method devised by Majorana 
leads to a semi-analytical series expansion, obeying both boundary conditions 
for a neutral atom 

0(0) = 1, (9a) 
0(oo) = 0. (9b) 



In a recent work, Guerra and Robotti jzzj have rediscovered a not well 
known short communication by Majorana, entitled Ricerca di un'espressione 
generate delle correzioni di Rydberg, valevole per atomi neutri o ionizzati 
positivamente ( Quest for a general expression of Rydberg corrections, valid 
for either neutral or positively ionized atoms) [7]. In that work, perhaps 
in the attempt of improving the asymptotic behaviour of the solution to 
Thomas-Fermi equation, Ettore requires that the potential vanishes for a 
certain finite value of x, say x , both for neutral atoms and for positive ions. 
He writes the self-consistent potential as 

ry 

V(r) = — (P + C, (10) 
r 

where, for an atom positively ionized n times (n = Z — N) , the constant C 
equals 



where 



1 f Z -n 



2/3 

» = 0M Mz^t) A ' (12) 



and the boundary conditions to Eq. (J3J) now read 



0(0) = 1, (13a) 
0(x o ) = 0, (13b) 

'(so) = (13c) 



One immediately notices that, due to the new boundary conditions, Eq. (fTU|) 
does not reduce to Eq. (jlj) even for n — 0, i.e. for a neutral atom. In other 
words, Majorana does not consider the potential V(r) in a generic location 
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Figure 1: Thomas- Fermi screening factor for the self-consistent potential of 
a neutral Ne atom (Z = N = 10). Solid line is Fermi's solution, dashed line 
is Majorana's solution, while the light dashed line has been obtained within 
Hartree-Fock approximation. 

of the electron cloud, but the effective potential acting on a single electron, 
thus excluding the interaction of an electron with itself. 

Probably, owing to his profound critical sense (let us remind that his 
colleagues in the Panisperna group nicknamed him the 'Great Inquisitor'), 
Majorana must have not excessively relied on his own solution [TJ2], which 
however reproduced the numerical solution of Thomas- Fermi equation quite 
accurately. Probably, Majorana was looking for a solution which should not 
decrease so slowly as x — > oo, as Eq. (JSJ) does. 

In Fig. ^ we report Thomas- Fermi screening factor as a function of r for 
a neutral Ne atom (Z = N = 10). The solid line refers to Fermi's numerical 
solution, with boundary conditions given by Eq(s). JHJ), the dashed line refers 
to Majorana's solution, with boundary conditions given by Eq(s). (|13|) with 
n — 0, while the light dashed line has been obtained within the Hartree-Fock 
approximation (see App.[X]for a derivation). As it can be seen, Majorana's 
solution introduces only a minor correction to Fermi's solution at finite x 
values, but is strictly zero for x > Xq. 

In his work on positive ions [23!) Fermi considers a potential vanishing 
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r[A] 

Figure 2: Thomas-Fermi screening factor for the self-consistent potential 
of the Ne + ion (Z = 10, iV = 9). Solid line is Fermi's solution, dashed line 
is Majorana's solution, while the light dashed line has been obtained within 
Hartree-Fock approximation. 

at a finite value x = Xq. However, instead of Eq(s). (fl"3|) . he employs the 
boundary conditions 

0(0) = 1, (14a) 
-x (j)'(xo) = — , (14b) 

which in particular imply Eq. (j9b|) in the case n = 0, corresponding to a 
neutral atom. 

In Fig. |2l we again report Thomas- Fermi screening factor as a function 
of r according to Fermi, Majorana, and Hartree-Fock, respectively, but now 
for a positively ionized Ne atom, Ne + {Z — 10, N — 9, n — 1). Majorana's 
solution again differs but marginally from Fermi's solution, but while for 
a neutral Ne atom Fermi's solution decreases too slowly, it decreases too 
rapidly for Ne + . 

Here, we are not disputing whether Majorana's note, Ref. [7], should be 
considered as a 'full' paper [21] , nor do we want to undervalue the importance 
of the contribution analyzed in Ref. ^H] ■ We would rather like to emphasize 
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that Majorana was conscious that his correction 1 did not lead to substantial 
modifications to Fermi's solution of Eq. ([5)1. including in the asymptotic limit 
(x» 1) j2Zj. 

Ettore never published anything else on this subject. 

4 Helium molecular ion 

In his successive work jc\, Majorana deals with the formation of the molec- 
ular ion He^. There again, Majorana demonstrates his exceptional ability 
to focus on the main physical aspects of the problem, while showing the 
limitations of his own theoretical approximations. He immediately observes 
that the problem is more similar to the formation of the molecular ion H^ , 
than to the reaction He + H. The most relevant forces, especially close to 
the equilibrium distance, are therefore the resonance forces, rather than the 
polarization ones. By exchanging the two nuclei, the system remains un- 
changed. Majorana makes then use of the method of Heitler and London 
[28] , and emphasizes the importance of inversion symmetry with respect to 
the middle point between the nuclei, set at a distance R apart. 

Heitler and London [2H] introduced a relatively simple expression for the 
wave-function \1/ of the two electrons in a hydrogen molecule H2 in terms of 
the wave-functions <p and ip of one electron in the atomic orbital correspond- 
ing to atom a and b, respectively: 

*(1,2) = (p(l)^(2) ±V(%(2), (15) 

where 1 and 2 denote the coordinates of the two electrons, respectively. The 
wave-function \I/s, corresponding to the choice of the plus sign in Eq. ([15)1. is 
symmetric with respect to the exchange of the coordinates of both electrons 
and nuclei, while (minus sign in Eq. ([15)1) is antisymmetric. The full 
wave-function is globally anti- symmetric, but here we are neglecting its spin 
part, since the Hamiltonian is spin independent. 

Despite its simplicity, the success of Heitler-London approximation relies 
on the fact that it explained the stability of the H2 molecule, and could 
reproduce with remarkable accuracy the dependence of the total electronic 
energy Ej on the internuclear distance R. One obtains the attractive solution 

1 Fliigge |2H] erroneously attributes this correction to Amaldi. Probably, he was only 
aware of Fermi and Amaldi's final work, Ref. |26| . 
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in correspondence with the eigenfunction \&s. It is relevant to stress at this 
point that, if one had considered only (p(l)ip(2), or ■0(1)^(2), in Eq. (fTo")). the 
agreement with experimental data would have been rather poor. Therefore, 
the resonance or exchange term is quite decisive for establishing the chemical 
bond. 

Heitler-London theory is even more accurate than the method of molecu- 
lar orbitals [2H1 EDI EU E2] (see, e.g., Ref. [33] for a more detailed discussion), 
which in addition to Eq. ()15|) takes into account also for the ionic-like con- 
figurations 

<p(l)<p(2) and 0(1)^(2), (16) 

corresponding to having both electrons on atom a, or b, respectively, on the 
same footing and with equal weights as the terms in Eq. (|15|) . However, 
the theory can be improved by adding to Eq. (|15|) the two contributions in 
Eq. ()16|) with appropriate weights, to be determined variationally. 

As in Heitler-London, in order to study the case of He^~, also Majorana 
starts from the asymptotic solution, namely for large values of R. In the case 
of H 2 , for large values of R, it is very unlikely that both electrons reside on the 
same nucleus. Similarly, in the case of He^~, Majorana neglects the possibility 
that all three electrons be located on the same nucleus. Ettore then proceeds 
by writing the unperturbed eigenf unctions for the three electrons (labeled 1, 
2, 3 below) in He^ as 

A 1 = ^ 23 , B 1 = 0!$ 23 , (17a) 

A 2 = ^ 2 *31, B 2 = 2 $31, (17b) 

A 3 = ^3^12, B 3 = 3 $ 12 , (17c) 

where $ and ip denote the wave-functions of the neutral and ionized a atom, 
respectively, while \1/ and ip denote the analogous wave-functions for atom b. 
Evidently, A 2 and A% can be obtained from A\ by permuting the electrons, 
and the i?'s from the A's by exchanging the nuclei. 

The interaction between the atoms mixes all these wave-functions, but 
by means of general symmetry considerations, first introduced by Hund j33] , 
as well as of inversion symmetry and of Pauli exclusion principle, Majorana 
concludes that the only acceptable wave-functions are 

Vl = A 1 -A 2 + B 1 -B 2 , (18a) 
y 2 = A 1 -A 2 -B 1 + B 2 , (18b) 
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which are antisymmetric in 1 and 2. In particular, y 2 corresponds to the 
(lscr) 2 2pa ( 2 E) configuration, viz. the bonding solution for the He^" molec- 
ular ion. The latter configuration is characterized by two electrons in the a 
molecular orbital built from the two Is atomic orbitals, one electron in the a 
molecular orbital built from the 2p atomic orbitals, as well as by a value of 
the total orbital angular momentum L = 0, and by a value of the total spin 
S = +|. The wave-function Eq. (j!8b|) clearly shows that the ground state 
is a resonance between the configurations He: — He- and He • —He: , where 
each dot denotes the presence of one electron on the a or b atom. 

In order to perform the calculation of the interaction terms explicitly, 
making use of analytic expressions, one can take the ground state of the 
helium atom as the product of two hydrogenoid wave-functions. However, 
it is well known that the result is greatly improved if, instead of taking 
the bare charge Z = 2 of the He nucleus, an effective nuclear charge Z e g 
is introduced, to be determined variationally. The fundamental effect here 
taken into account by Majorana is that of screening: In an atom with many 
electrons, each electron sees the nuclear charge Ze as slightly attenuated by 
the presence of the remaining electrons. 

The concept of an effective nuclear charge, already introduced for the 
helium atom, had been extended by Wang 35 to the hydrogen molecule. 
Probably Majorana was not aware of Wang's work, since he does not refer 
to it in his 1931 paper. In any case, Majorana is the first one to make use of 
such a method for He^J\ In making reference to his own work [HE], where Z e g 
is used as a variational parameter for He^", Pauling reports in a footnote 2 

"The same calculation with Z e g given the fixed value 1.8 was 
made by E. Majorana jS]." 

The variational value obtained by Pauling for Z e g is 1.833. 

By making use of his results, Majorana evaluates the equilibrium internu- 
clear distance as d — 1.16 A, in good agreement with the experimental value 
1.087 A. He can then estimate the vibrational frequency as n = 1610 cm , 
which he compares with the experimental value 1628 cm -1 . Majorana con- 
cludes his paper by stating [H] that his own result is 

"casually in perfect agreement with the experimentally deter- 
mined value" 

2 See footnote on p. 359 of Ref. EH]- 
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Figure 3: Variational energies of the molecular ion He^~, as a function of 
the internuclear distance R. Solid line refers to the symmetric wave-function 
in Eq. (|T3j). dashed-dotted line to the antisymmetric one, while dashed line 
refers to the 'non-bonding' case, where position exchange is neglected. Re- 
drawn after Ref. [HE] (see also Ref. 



(our italics). Any other author would have emphasized such a striking agree- 
ment as a success of his own method, whereas Majorana rather underlines 
the drawbacks of his own approximations. 

We would like to remind that he also estimates the minimum energy, i. e. 
the dissociation energy, finding the value E min = —1.41 eV, but he had no 
available experimental data to compare with, at that time. However, he is 
not satisfied with such a result and collects |H] 

"all the errors of the method under the words ' polarization forces'," 

which he estimates for very distant nuclei using the polarizability of the 
neutral He atom. He then finds E m[n = —2.4 eV. More recent theoretical 
calculations, using the method of configuration interactions [39J or ab initio 
variational methods |2Qj, have estimated the value E min = —2.47 eV. The 
experimental value has been accurately determined quite recently |H] as 
-Emin = -2.4457 ± 0.0002 eV. We are not claiming that Ettore's result is 
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more accurate than the theoretical results mentioned above. However, he 
certainly understood the essential physical effects for that system, and made 
use of appropriate approximations to estimate them. In particular, it is 
interesting how he emphasizes the quest for the symmetries of the system (see 
the translation of a paper by Majorana in Ref. |20j). As in the case of H2, also 
for He^~ it is essential to include the position exchange term between He and 
He + , in order to have chemical bonding, as it can be seen in Fig- El redrawn 
after Ref. [3H] . If one had neglected the resonance He : He + # He + : He (see 
dashed line in Fig. EJ), chemical bonding would have been impossible. 

5 The discovery of the neutron 

Rutherford's pioneering work [4j| paved the way not only to Bohr's atomic 
model, but also to nuclear physics. 

In 1930 Bothe and Becker [33], like Rutherford, employed a particles 
against a berillium target in a scattering experiment. They observed the 
emission of a very penetrating radiation, which they interpreted as 7 rays. In 
successive experiments, Irene Curie and Frederic Joliot |33] , her husband, 
developed further these experiments, but they arrived at similar conclusions. 
According to Emilio Segre's account 0H]> Majorana thus commented the 
Joliots' results: 

"They haven't realized they have discovered the neutral proton." 

At this point we should remind that at that time it was believed that 
the nucleus was composed by protons and electrons. It was Chadwick [37] 
who soon after demonstrated that the radiation emitted in the Joliots' ex- 
periments was made up by neutral particles, whose mass is very close to the 
proton's mass. It was probably Fermi [H)] who first distinguished between the 
neutrinos conjectured by Pauli, and the neutrons discovered by Chadwick. 

Meanwhile, Majorana developed a theory of the nucleus containing pro- 
tons and neutrons and then, according to Segre [35] . 

"he analyzed, as far as it was possible, the nuclear forces on the 
basis of the available experimental results, and he estimated the 
binding energies of the lightest nuclei. When he presented his 
work to Fermi and ourselves, we immediately recognized its im- 
portance. Fermi encouraged Majorana to publish his own results, 
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p n " p 

Figure 4: Exchange interactions. Resonant forms in the hydrogen molecular 
ion, H: H + ^ H + : H (upper row), and in the proton- nucleon pair inside a 
nucleus, p: n # n: p (lower row). 



but Majorana refused to do so, saying they were yet too incom- 
plete." 

More than that, when Fermi asked Majorana whether he could make refer- 
ence to his results during a forthcoming conference in Paris, Ettore mockingly 
replied he would agree, provided the reference was attributed to an old pro- 
fessor of electrochemistry, who was also going to attend the same conference. 
Obviously, Fermi could not accept Majorana's condition, and no reference 
was then made to his results during the conference. 

Meanwhile, people started feeling the lack of a theory of nuclear forces, 
conveniently taking into account for the presence of both protons and nucle- 
ons in the nucleus. But where to begin with? 



6 Heisenberg-Majorana forces 

To this aim, in three fundamental contributions |2HJ ESI EH] > Heisenberg as- 
sumed hydrogen molecular ion as a model. He recognizes that the most 
important nuclear forces are not the polarization forces among the neutrons, 
or Coulombic repulsion among protons, but the exchange forces between 
protons and neutrons. 
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Figure 5: Kinetic, potential, and total energies for the ground state of H^, 
excluding nuclear repulsion, within the linear combination of atomic orbitals 
(LCAO) approximation. Cf. Fig. 2.4 in Ref. [ST], where the same quantities 
have been obtained within a variational method. 

Heisenberg emphasizes that neutrons obey to Fermi statistics. Moreover, 
since a neutron possesses spin it cannot be simply thought of as composed 
of a proton plus an electron, unless the latter has zero spin, when inside a 
neutron. 3 A neutron is an elementary particle per se. The interactions 
postulated by Heisenberg are characterized by the exchange of both position 
coordinates and spins of the two nucleons. 

Similarly, Majorana assumed that the fundamental nuclear forces are of 
exchange nature between protons and neutrons. However, he fully exploits 
the analogy with H^ (see Fig. HJ, regardless of spin. 4 

3 Besides considerations concerning the spin, such a model would require an enormous 
amount of energy to localize the electron within the neutron |49j . 

4 Current literature usually employs the formalism of isotopic spin to describe the ex- 
change character of the nuclear forces. However, as noted by Blatt and Weisskopf |52j. 
this is equivalent to a description which makes use of the forces of Bartlctt, Heisenberg, 
Majorana, and Wigner. 
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Let a i and r 2 , cr 2 stand for the position and spin coordinates of the 
first and the second nucleon, respectively, and let i/j(r%, o\\ r 2 , a 2 ) be the wave- 
function for a given nucleon pair [521. Then Heisenberg exchange P H implies 

P H tp{r 1 , <7i; r 2 , cr 2 ) = tp{r 2 , cr 2 ; ri, 01), (19) 
whereas Majorana exchange P M implies 

P M ^(ri, <ri; r 2 , a 2 ) = ^(r 2 , m; r 1; a 2 ). (20) 

In Majorana's own notation (apart from a minus sign here included in the 
definition of J(r)), the exchange interaction then reads [2] 

(Q', q'\J\Q", q") = J(r)6(q' - Q")5{q" - Q'), (21) 

where Q and q are the position coordinates of the neutron and the proton, 
respectively, and r = \q' — Q'\ is their relative distance. Majorana then plots 
a qualitative sketch of J(r) (cf. Fig. 2 in Ref. [2]), which closely resembles 
the behaviour of the potential energy in fhj" , when the internuclear repulsion 
is neglected (Fig. EJ). 

In the same paper |2] , in addition to his knowledge of molecular physics, 
Majorana fully exploits also his acquaintance with the atomic statistical 
model. Indeed, he defines the nuclear density as 

8tt 

w 

in complete analogy with Eq. (JTJ), where P n and P p are the Fermi momenta 
of neutrons and protons, respectively. From this model, he derives an asymp- 
totic expression (p — ► oo) for the exchange energy per particle, 

o( " )l -» = -^ J(0) ' (23) 

where n\ and n 2 are the numbers of neutrons and protons, respectively. As 
in Thomas-Fermi model, the kinetic energy per particle, t say, is given by 

t oc p 2 '\ (24) 

From the competition between kinetic and potential energy, the total energy 
attains a minimum as a function of r (cf. Fig. 1 in Ref. [2]). 

Majorana's model explains two fundamental properties of nuclear physics 
[221 : (a) the density of nucleons is about the same for all nuclei {density 
saturation)] (b) the binding energy per nucleon is about the same for all 
nuclei {binding energy saturation). 



P=^( p n + P*l (22) 
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7 Concluding remarks 



In some of his fundamental papers, Majorana mainly focussed on the asymp- 
totic properties of the potential and of the wave-function of an atomic or 
molecular system. This is clearly demonstrated in his work on helium molec- 
ular ion, He^~ jB]. On the basis of his hypercritical spirit, Majorana was 
probably unsatisfied with the asymptotic behaviour of the screening factor 
(f) within Thomas- Fermi model, but his note [7j is too short to confirm that. 
What we can certainly emphasize is his taste for the quest of symmetries, 
and their relevance to determine the main properties of a physical system 
[21]. This led him to demonstrate that the exchange symmetry is essential 
to the formation of the chemical bond. Exchange symmetry is also central 
in his model of the nuclear forces. 

The quest for symmetries is evident in his famous work on the symmetrical 
theory of the electron and the positron j^j. There, he notes that 

"all devices suggested to endow the theory with a symmetric 
formulation, without violating its contents, are not completely 
satisfactory. [. . . ] It can be demonstrated that a quantum the- 
ory of the electron and the positron can be formally symmetrized 
completely by means of a new quantization procedure. [. . . ] Such 
a procedure not only endows the theory with a fully symmetric 
formulation, but also allows one to construct a substantially new 
theory for chargeless [elementary] particles (neutrons and hypo- 
thetical neutrinos)." 

Several important experiments jSJ E] are currently under way to observe 
the 'Majorana neutrino'. 
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A Thomas- Fermi screening factor 

within Hartree-Fock approximation 



In order to critically assess the accuracy of Fermi's and Majorana's approx- 
imate solutions for the atomic screening factor <p in Thomas-Fermi model, 
let us briefly derive it within the Hartree-Fock self-consistent approxima- 
tion. The solution of Hartree-Fock equations enables one to determine the 
(spherically symmetric) radial electron density 

D(r) = Anr 2 p(r), (25) 

normalized to the total electron number as 

POO 

/ D(r) dr = N (26) 
Jo 

(see, e.g., Ref. [HI]). The radial electron density D{r) of a neutral Ne atom is 
characterized by two peaks, referring to the Is and 2s 2p shells, respectively 
(cf. Fig. 8.6 in Ref. jM])- In the case of Ne + , such peaks are slightly shifted 
at smaller values of r. Although D(r) is always strictly different from zero 
over the whole r range, it is an exponentially decreasing function of r, with 
D(r) ~ roughly defining the atomic (respectively, ionic) radius. 

By relating the electric field \E\ = (l/e)dV/dr corresponding to the self- 
consistent potential, Eq. (jlj), to that generated by the nucleus and the elec- 
tron cloud within a distance r from the nucleus, by Gauss law, one finds 



0(0) 



1 
r 
1. 



1 + 4 / D(r')dr' 
" Jo 



(27a) 
(27b) 



where a prime here refers to derivation with respect to r. 

Within Thomas-Fermi approximation, 0(r o ) = 0, where r is the ionic 
radius, and the integration in the normalization condition, Eq. (j26j) . should 
actually be performed up to r = Tq. Then, Eq. (127a)) reduces to Fermi's 
boundary condition, Eq. ()14b)) . 

Within Hartree-Fock approximation, D(r) is in general nonzero at any 
finite r. However, as r — > oo, the potential experienced by a test charge is 
that of a charge (Z — N)e, i.e. V(r) — Ep ~ —(Z — N)e 2 jr. Comparing such 
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an asymptotic behaviour of the potential with the definition of the screening 
function in Eq. (@J), one has 

N 

lim 0(r) = 1 - — . (28) 



On the other hand, making use of the latter result, from Eq. ()27aj) it follows 
that 

lim r0'(r) = 0, (29) 

r— >oo 

thus implying that <p'(r) vanishes as r — > oo more rapidly that 1/r (in fact, 
it vanishes exponentially). 

Finally, from Poisson equation, — V 2 V = 4ire 2 [—Z5(r) + p(r)], for a given 
electron charge distribution, Eq. (}25|) . one obtains 

<p = --U. 30 
z r 

at any r > 0. Integrating once between r and oo, and making use of the 
limiting value Eq. ()29j) . one obtains 

1 r°° D(r') 

4f(r) = -± ^V-dr', (31) 



Z . 

which combined with Eq. (|27ajl yields the desired screening factor 

1 r r r°° D(r') 

0(r) = 1 - - ^ D(r') dr' - - y dr' (32) 

in terms of the Hartree-Fock self-consistent radial density D(r). In particular, 
Eq. (j%2} manifestly fulfills the boundary conditions 

0(0) = 1 (33a) 
N 

0(oo) = 1--. (33b) 

In Fig(s).dand|H light dashed lines represent Eq. (|32|). with D(r) numerically 
obtained within Hartree-Fock self-consistent approximation for Ne and Ne + , 
respectively. 
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